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We study the liquid-glass transition of the Lennard-Jones binary mix- 
ture introduced by Kob and Andersen from a thermodynamic point of view. 
By means of the replica approach, translating the problem in the study of a 
molecular liquid, we study the phase transition due to the entropy crisis and we 
find that the Kauzmann's temperature Tk is ~ 0.32. At the end we compare 
analytical predictions with numerical results. 
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I. INTRODUCTION 



paper; 



The aim of this paper is to study the liquid-glass transition and investigate the glass phase of 
a binary mixture of particles which describes a fragile liquid. We consider a binary mixture 
interacting via a Lennard- Jones potential, introduced by Kob and Andersen0. This model is 
suitable for such problems because of the lack of a crystalline state. This allows a detailed 
comparison between analytical results and numerical simulations. 

The proposed theory relies on the following assumptions, discussed in detail in previous 
sH: 

In the supercooled liquid phase, the phase space can be partitioned in an exponentially 
high number of potential energy minima separated by energy barriers of order 0(1) (the 
barriers will diverge at a low temperature, at the Kauzmann temperature where the 
viscosity diverges). The basins of attractions of these minima are the so-called 'inherent 
structures'!. A growing number of numerical studies are supporting this picture@~0. 
Choosing the thermodynamic point of view, we focus mainly on free-energy landscape. 
We assume that below the dynamical (Mode CouplingS) critical temperature, T^BEl, 
the supercooled liquid at equilibrium is almost always trapped in one of the exponentially 
large number of local free energy minima, which we will call 'valleys'. This picture is 
very similar to the one which describes generalized mean field models of spin-glasses 
displaying a one step replica symmetry breaking transition0^0. The free-energy minima 
in short range systems are not separated by infinite barriers, so different valleys with 
the same free energy can be explored by the system during its evolution. This is what 
is supposed to happen in the liquid phase. 

We write the partition function as a sum of contributions from different minima: 



e -AW>(T) 



/ rf/e-^A/m (1) 



where J\f(f, T) is the number of minima of the free energy at temperature T as function 
of the free energy density /. We can show that the system does not choose the con- 
figuration with minimal free-energy, but instead minimizes the generalized free energy 
density, namely: 

0(T)~min[/-TS(/,T)]. (2) 



The total entropy density is given by the sum of the entropy density of a typical minimum 
and of the complexity, or configurational entropy £ ~ logAf/N, taking into account the 
number M of these minima. 

• In the glass phase, the system is described, to a first approximation, as an amorphous 
solid, where the only degrees of freedom are small displacements around disordered cen- 
ters of oscillation. We assume further that thermodynamic quantities are self-averaging 
with respect to the disorder. 
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The main hypothesis relies on old Adam-Gibbs-Di Marzio scenario where the liquid- 
glass thermodynamic transition is related to an 'entropy crisis', that is the existence of a finite 
temperature where ergodicity breaks and the configurational entropy vanishes. Within this 
picture, the transition temperature, TV, corresponds to the temperature where the viscosity 
diverges with a generalized Vogel-Fulcher law!, 7] oc exp(T — Tk) v ■ In other words, the typ- 
ical singular behaviour of experimental measures of thermodynamic quantities (e.g. specific 
heat) at a somewhat arbitrary T g is related at a 'true' second order transition which occurs 
at a lower temperature. The internal energy is continuous at Tk, there is no latent heat and 
the specific heat suddenly falls from the liquid value to a definitely smaller one, in agreement 
with the Dulong and Petit law. Nevertheless, from the point of view of an appropriate order 
parameter, we find that the transition is discontinuous. Indeed, while in the glass phase each 
particle is confined in a 'cage' of finite size due to interactions with its neighbors, in the liquid, 
because of diffusion, there is no cage, i.e. there is no thermodynamically stable cage, but only 
some metastable 'caging' effect. 

The glass phase is studied by generalizing a recently proposed theoryi'S to binary mixture, 
allowing one to obtain equilibrium properties of fragile glassesil from those of liquid phase, 
computed for a molecular liquid consisting of m 'clones' of the system with! m < 1. This 
extension has been accomplished so far only for a soft sphere binary mixture mo del, while 
in the present paper we address the more realistic model of a Lennard- Jones binary mixture. 
Due to the different nature of the system and to the presence of a liquid-gas transition, the 
simple Hyper-Netted-ChainfH approximation previously used! is not enough precise. We must 
introduce some slightly different approximation in investigating liquid and glass phases which 
will be discussed in detail later. 

After first introducing the model (section II), we present and discuss approximations 
needed to compute replicated free-energy (sections III-IV), we then describe a way to study 
supercooled liquid phase (section V) and finally we present analytical results and comparisons 
with numerical simulations (VI). 



II. LENNARD-JONES BINARY MIXTURES 

We study mixtures of two types of particles with different radii, here called + and — , with 
pairwise interactions. The Hamiltonian of our problem is: 

H= £ I'" (.r,-,- ; ) e, €{-,+}, (3) 

l<i<j<N 

where the particles move in a volume V of a <i-dimensional space, and V ++ , V + ~ , and V 
are arbitrary short range potentials. 

The so-called soft sphere (SS) model, where the potentials V(r) are short ranged and 
purely repulsive (i.e. proportional to R~ 12 ), was studied in a previous work!. In that paper 
one: 

• demonstrated the feasibility of extending the thermodynamic theory of liquid-glass tran- 
sition recently presented™ to binary mixtures; 

• compared analytical results to numerical simulations; 
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• compared different methods of approximation proposecffl to compute 'molecular' or 
replicated free-energy. 

The free energy in the uncorrelated liquid phase and in the correlated glass phase were 
there computedi using the framework of HNC closured. By the HNC approximation the 
computed entropy in liquid phase was found correct to within 10 %. Other approximations 
were introduced to compute the entropy in glass phase, namely the 'quenched' approximation, 
which neglects the feedback of vibrational modes on center of masses, and the 'superposition' 
approximation, where p-point correlation functions are expressed as 'chain' products of two- 
point correlation functions. When more realistic potentials are introduced, it is sometimes 
necessary to use more sophisticated integral equations than those based on the HNC. 

In this paper we study a model of binary mixture introduced by Kob and Andersen0, 
where potentials are of the Lennard- Jones (LJ) type: 



V €e (r) = 4rf 



a ff < \ 12 / o ff , 



r 



(4) 



The concentrations used are c+ = | and c_ = |. With the choice of parameter values 

a ++ = 1.0, cr + _ = 0.8, a = 0.88, r] ++ = 1.0, r/ + _ = 1.5 and r\ = 0.5 it is possible to 

prevent both crystallization and separation in two phases with particles of different type. For 
the sake of clarity we underline that the symbols + and — correspond to A and B of the Kob 
and Andersen paper and, more importantly, here we do not use any cut-off on the potentials. 

The main effect of the introduction of an attractive term in potential V(r) is the appear- 
ance of a density-dependent temperature, Tql{p), where the system undergoes a separation 
into a liquid and gaseous phases with different densities. 

This model has been extensively studied numerically in the past0&~§II>EI. In particular, 
at density p = 1.2, it has been shown that a 'dynamical' transition occurs at temperature 
Td = 0.435. Below this temperature aging phenomenaBHEl] and violation of the fluctuation- 
dissipation theoremcZH are observed. 

In this paper, the same model is studied from a 'thermodynamic' point of view, showing 
the existence of a temperature, T K = 0.32 (< T D ), where the system undergoes a real liquid- 
glass transition. 

As a preliminary step, let us note that we have to go beyond HNC approximation to 
study the model in the fluid phase. Indeed, within this closure, integral equations at density 
p = 1.2 show a singularity at T GL ~ 0.5, overestimating the true T GL . We are forced, 
then, to introduce some more reliable approximation in order to study this model. The 
correlation functions used in this paper are computed solving an integral equation which 
continuously interpolates between the HNC approximation at long distances and the Mean 
Spherical Approximation (MSA) at short distances, as proposed by Zerah and Hansenlll. The 
details are described below. 



III. GENERAL FRAMEWORK 

The basic ideas used to study the glassy phase are the ones discussed in the introduction 
and in previous papersi~i. We briefly review here how to use them for explicit computations. 
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Letting f eq be the equilibrium free energy, which minimizes 0, the system enters the glass 
phase when the complexity S(/ eg , T) vanishes. Therefore the computation of £(/, T) it is the 
first fundamental topic of our investigation. 

This issue can be settled by using the method proposed by MonassonS, where m identi- 
cal replicas of the original system are considered. By introducing a small but extensive and 
attractive coupling term among them, these replicas are forced to be in the same valley, so 
that they become strongly correlated. In the limiting case of no coupling term, two arbitrary 
replicas in the liquid phase are always found in two different valleys, because of the exponen- 
tially large number of them, which makes the probability of finding any two systems in the 
same valley negligible. The transition to glass phase is supposed to be related to an ergodicity 
breaking, where the free-energy barriers become infinitely large. In the glass phase replicas 
remain correlated even after sending the coupling term to zero. In a sense, the coupling term 
acts like, for example, the external magnetic field in the Ising model, exposing the transition. 
Moreover, the crucial point is that using the same argument as for (^|) the free energy of the 
replicated system <3>(m, T) is related to the complexity £(/, T) by the simple formula: 

$(m,T) ~nun[m/(T) -TE(/,T)], (5) 

It is clear from this formula that if we are able to compute this thermodynamic potential 
for an arbitrary, non-integer value of the parameter m, we can obtain all the values available 
for free energy, /, and complexity, S, simply by varying this parameter m in the following 
formulas: 

d$(m,T) _ j m 2 d($(m,T)/m) _^ 
dm T dm 

By eliminating m from these formulas one obtains the complexity as a function of the free 
energy £(/). Obviously, the physical free-energy and configurational entropy are obtained in 
the limit m — > 1. 

Because the function $(m, T)/m should be a convex function of parameter m, a physically 
possible situation is described only in the interval m e [0, m*], where m* is the point of 
maximum. Indeed, m* corresponds to the minimum value f m in available for free-energy. For 
greater values of m, the second equation in implies that complexity is negative, which 
cannot be accepted on physical grounds. 

The region of temperatures where m* > 1 is the liquid phase, where equilibrium complexity 
is positive and the free energy is that which minimizes the balance between f3f and £(/, T). 
In the low temperature glass phase, m* < 1 instead, and in the whole region m e [m*,l] 
the replicated free energy (per replica) has the value $(m, T)/m = f m in(T). In particular we 
have that the equilibrium free energy $(1,T) = f m in{T). From this free-energy it is possible 
to get all the other thermodynamic quantities, internal energy, specific heat and so on in the 
usual ways. 

It is useful to think about replicated free-energy as the free energy of a molecular liquid, 
whose molecules are composed of m atoms, each in a different replica of the liquid. The 
tendency to form molecules is forced by a small but extensive coupling term between particles 
of different replicas@'E3. 

In the case of binary mixtures, one may ask if these molecules can be formed by particles 
of different kind. For the soft sphere binary mixture mo dell the point is quite subtle because 
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the particles (and interactions between them) are different only because of their effective size, 
and the ratio of their radii is not so different from 1. In the Lennard- Jones binary mixture 
model the situation is more clear. Beside the sizes of particles, there are three different values 

of the potential well for the three different kinds of interaction (++, , H — ). In this case, 

the energy depends much more strongly on the type of particles involved. In fact, we expect 
that even the substitution of a single particle of a given type with a particle of different 
type, at low temperature, substantially increases the total energy of the system, lowering the 
probability of such a configuration which can then be safely neglected. 

Since replicated (or molecular) free energy must describe the free energy of m systems in 
the same state, we introduce attractive coupling terms only between particles of the same kind, 
neglecting molecules formed by atoms of different types. The replicated partition function is 
then: 




- E E w +(^«-^('))- E - x * b (i)) ] > ( 7 ) 

ig{+} a^b «G{-} a + h j 

where uo is the attractive potential among replicas and we must sum over the permutation (a a 
and 7r a ) of the particles in each replica. When relabeling particles, so that to particle i of a 
given kind in replica a corresponds particle i of the same kind in replica b (which is supposed 
to belong to the same molecules) and so on, the sum over permutations can be eliminated 
leaving a multiplicative factor (N + \ AM)^ -1 ). 

The_phase transition we are investigating is equivalent to that shown by generalized spin 
glasses^. The main idea is that in the glass phase there are only few states, not linked by 
trivial symmetry transformations, which contribute to integral ([!]). Below Tk these states 
are mutually inaccessible because of ergodicity breaking. In replica space, this spontaneous 
one-step replica symmetry breaking is signaled by the onset of an off-diagonal non trivial 
correlation, when couplings are sent to zero. 

The study of the thermodynamic phase, where replicas are correlated, is accomplished by 
choosing the probability distribution (after selecting a permutation) of the m particles inside 
the same molecule for both kinds of particles as order parameters: 

p + (r\...,r m )= £ <8{x\-r 1 )...8{x™-r m )> 
ie{+} 

p_(r\...,r m )= £ <5(x 1 l -r 1 )...5(xT-r m )> (8) 
ie{-} 

where a permutation is already selected, and introducing the Legendre trasformate of molec- 
ular (replicated) free energy: 

G m [ P }= lim l-J-logZ m [u;]-- [f[d d r a £ p £ (r\ ...,r m )W e {r\ ...,r m ) ) (9) 

,.. _^ n \ pm mJ a=1 e =+- / 
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with 



W e (r\...,r m ) = Y,^(r a -r b ) (10) 

a<b 

Performing the limit uj — > is equivalent to search for a non-trivial saddle point, with respect 
to the order parameters, of the functional G m [p]. In the presence of a glass transition we 
expect the order parameters and thermodynamic quantities to show the following behaviours: 

• For T > Tk the density free energy is the liquid one (m = 1) and the order parameters 
are trivial, i.e. p + = c + p and p_ = c_p (where the c's are the different concentrations). 

• For T < Tk, the density free energy of glass is the maximum of function f m = 
$(m, T)/m, that is found at m* < 1. From f m * we can compute all the thermody- 
namic quantities. 

The free energy and its first derivatives are continuous at TV, while the heat capacity 
falls suddenly to solid-like values. The transition, then, is of second order from thermo- 
dynamics point of view, but it is discontinuous in the order parameter which abruptly 
becomes a non trivial function of positions in different replicas. 

It is natural to describe the particle positions in term of center of mass coordinates and 
relative displacements uf with x\ = Zi + uf and J2 a — 0. A useful simplification is the 
choice of quadratic coupling among replicas, that allows rewriting (0) as: 

N / m \- 

n U d *(£«?) 

,i=l \ a=l / . 

• exp ( -13 E V ^ ( z i ~ z i + M i ) 

\ a=l i<j 

-^EEK-^-^EEK-^) 2 ), (") 

a,b ie+ H:LX ~ a,b ie- / 

where the {u^} for a given i (the indices p and u, running from 1 to d, denote space directions) 
turn out to be Gaussian random variables with a vanishing first moment and a second moment 
given by 

KO= W«— • (12) 

IV. REPLICATED FREE-ENERGY 

Three different ways to compute the 'molecular' free energy of a replicated system were 
introduced^. In a previous paper! on binary glasses, free energy was computed by two of 
these schemes of approximation, namely harmonic resummation and small cage approxima- 
tion; both gave similar results. Here we want to unify these two approximations in a single 
framework. 



7 



In the glass phase molecules are expected to have a small radius, so our starting point is a 
quadratic expansion of V in the partition function ([11]). The integration over these quadratic 
fluctuations gives: 

Nd/2 /n~ Nd ( m - 1 ) .TV / 1 



Z m = Ar ;, Ar _, / II ^ exp -/?m 2 ^ - - -^-Tr log (/3M) (13) 



i=i \ «<j 

where the matrix M, of dimension Nd x iVcf, is given by: 



m 



v z ( - z fc ) + — - V*? { Zi - Zj ) (14) 

e i / 

and V^ v (r) = d 2 V/dr fl dr 1/ . We have then an effective Hamiltonian where the centers of 
masses Z{ of the molecules interact at effective temperature T* = l/(/3m) by means of a 
pair potential, complicated by the contribution of vibration modes. We proceed by using a 
'quenched approximation', i.e. neglecting the feedback of vibration modes onto the centers of 
masses: 



Z m = m Nd l 2 V^ Nd{m ~ l) Z Hq ^m) (exp (~^^ \og^M{z) 



(3 m 



~ m Nd /^ Nd(m ~ 1) Z hq ((3m)exp (Tr log/3M(z)), m ) . (15) 

This approximation becomes exact when close to the Kauzmann temperature where m — > 1. 
We normalize matrix elements as follows: 



s y if ^ (*-%)• ( l6 ) 

introducing the mean values of diagonal terms: 

r e = J2 CeP I d d rg e€ ,{r)-AV~\r) + -. (17) 
The replicated free energy is: 

A M m — 1) 1 



(/ ' 1 ^ (c+ log(/3 r+) + c_ log(/3 r_)) 



2m 
1 (m-1) / 



(18) 



We consider fluctuations of diagonal term J2kCik up to second order, while the whole non- 
diagonal contribution is resummed by means of a kind of chain approximation. The 
free-energy to compute is, then: 
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(3) = ~ log(m) - ^L-il log( 27 r) - -L log Z Hq ((3 m) 
+ d {m 2 ~ l) (c + log(/3 r + ) + c_ log(/3 r_)) 



1 (m-1) 
~ 2N 2m 

The p-th order term depends, as usual, on the p-point correlation function 

(Trcn = £ E / • • • d " z p (fg^izi ■ ■ ■ z P )c;\%( Zl - z 2 ) 

ei...e p €{+,-}Mi-^P 

• • -C%-_\U z v-i ~ z P )C%U z p ~ (2°) 

The 'chain' approximation in the computation of the traces neglects contributions coming from 
the terms where two indices coincide and utilizes the superposition approximation to obtain 
p-point correlation functions g^ p \zi...z p ) = g{z,\ — z 2 ) • • ■ g(z p — zi). With these hypotheses 
we arrive at: 

(TrC) = ( d d Zl ...d d z p( f ]T E 9 £1£2 (z 1 -z 2 )C;\%(z 1 -z 2 ) 

Hi..-Hp ei...e p 

• • -g^{z p ^ - z p )C%-_\X(z P -i - z p )g^(z p - z 1 )C;% 1 (z p - z x ). (21) 

Traces are computed in Fourier space, where the matrix is diagonal with respect to particle 
indices, by introducing 

D%{k) = J d d rg~'(r)C£(r)e ikr , (22) 

which can be decomposed in its diagonal (longitudinal) and traceless (transversal) parts with 
respect to spatial indices 

D%{k) = ^ a-'(k) + {^- 6 -f) b-'(k). (23) 

The last step consists in the diagonalization of D in the space of components. For each k, 
one has four distinct eigenvalues, the two 'longitudinal' ones, corresponding to those of 

Df'(k) = a te \k) + ^r^b ee '(k), (24) 
and the two 'transverse' eigenvalues of the matrix 



Df'(k)=a ee '(k)-h ee '(k). (25) 



The eigenvalues are: 
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1 
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\± = 
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1 
2 




+ D 



+ v^ii 
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— 


\2 i A ( rH — \2\ 

) + 4 l- L 'll J J 


- y/iPV 




) 2 + 4(D+-) 2 ) 


+ v^i" 




) 2 + 4pr) 2 ) 






) 2 + 4( J D+-) 2 ) 



(26) 

The higher order terms involving the fluctuation of J2k Cik cannot be resummed in a systematic 
way similar to the one introduced for the off-diagonal terms . In any event, neglecting them 
is not likely to dramatically change the final result, because the sum of all terms higher than 
order two contributes just a small correction to second order computation for non-diagonal 
matrix elements. This is confirmed by comparison to simulations. 

In addition, the second order term from these calculations is equal to the second order 
term computed in the small cage approximation when non-harmonic corrections due to higher 
derivative are neglected. Summarizing: 

• The harmonic resummation schemeffl! allows us to resum the whole contribution com- 
ing from off-diagonal matrix elements, but it misses contributions from non harmonic 
corrections and from thermal fluctuation of diagonal elements of M. 

• Both kinds of correction are accounted by the small cage approximation!!!, but in that 
case we are unable to systematically resum either the terms of the series: 

tl^Y (27) 



or the contributions from non-harmonic corrections, and one has to truncate the expan- 
sion at some order (up to now we truncated at the second order). 

• It is possible to partially merge second order contributions from the small cage compu- 
tation in the harmonic resummation scheme, obtaining the whole resummed series plus 
fluctuations of diagonal elements up to second order. 

The final expression we get for the binary mixture free energy within this scheme is: 

// n\ d , i \ d(m — 1) . , n . dim — 1) , . ,„ 

<f>{m, (3) = -— log m - -i- >- log(27r) + ^— >- c+ log(/3r+) + c_ log(£r_) 

2m 2m 2m 

- ^rl / £ ."WWW + (i - i) *WW<*) 

r ee'e" 

{m ~ l) J d d rpY J 9 e "{r)Y J {C%{r)f - -L log Z liq (f3 m) , (28) 



2m 
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where the function L3 is log(l — x) + x + x 2 /2. Let us point out that the term in the second 
line is just the contribution, to lowest order, due to fluctuations of the diagonal elements of 
the Hessian matrix. 

It is now simple to compute the equilibrium complexity E eg : 



E eq (/3) = m 



2 d(3F n 



dm 



— Sli q — Sgol, (29) 

m=l 



where Su q is the entropy of the liquid at the effective temperature T e ff, which equals T for 
m = 1, and S so i is the entropy of a harmonic solid with a matrix of second derivatives given 
by M, i.e: 

S sol = d - log(27r) - ^(Tr log (/?M)>. (30) 

The condition for identifying the Kauzmann temperature, E eg (/3) = 0, reads in our ap- 
proximation simply: 

Sliq = S so l, (31) 

as expected on general grounds!!. 

If from our computation we find that Su q < S so i, we are already in the glass phase (T < Tk) 
and the computation in the liquid does not make sense, while in the other case Su q > S so i, the 
temperature is greater than T K (and of course less than T D if the spectrum of M is positive). 

With Q2"8|) we have a tool to investigate the thermodynamics of the low temperature glass 
phase. It should be emphasized that only the liquid phase g(r) and free-energy are needed 
in order to compute the glass phase thermodynamics. Beyond the usual thermodynamic 
quantities (energy, entropy, heat capacity...), we are interested in the two new parameters 
describing the glass phase: 

• The square cage radii A e , defined as A e = |((a^) — (xi) 2 ) for type e particles. These 
square cage radii are obtained by differentiating the free energy with respect to coupling 
terms and by sending couplings to zero at the end: 

A * = Tt 2 iw WTTste = (32) 
d(m — l)N e 0(1/ a e ) 

This expression gives square cage radii nearly linear in temperature in the whole glass 
phase (see figure [6b]), which is natural since non-harmonic effects have been neglected. 
The values of A e at the Kauzmann temperature are of the same order of those obtained 
for the soft sphere model, i.e. ~ 10 -3 . 

• The effective temperature T e // — T/m* of the molecular liquid. This varies very little 
and remains close to the Kauzmann temperature throughout the low temperature phase, 
confirming the hypothesis that glass can be successfully described by means of a system 
of molecules remaining in the liquid phase. We stress that the linear behaviour of 
parameter m* as a function of T is a feature shared by every 1RSB system to our 
knowledge. This static parameter should be equivalent to the dynamical parameter 
that measures violation of the fluctuation-dissipation theoreml!, as it is observed in 
short range spin-glasses0 
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The harmonic expansion makes sense only if M has no negative eigenvalues. This is 
natural since it is intimately related to the vibration modes of the glass. Notice that we 
cannot describe activated processes here, and so we cannot see the tail of negative eigenvalues 
(with number decreasing as exp(— 1/T) at low temperatures), which is always present. It is 
known however that the fraction of negative eigenvalues of M becomes negligible below the 
dynamical transition temperature, TdSEI. Therefore our harmonic expansion makes sense if 
the effective temperature, T e ff, is less than Tp. 



V. INTEGRAL EQUATIONS FOR BINARY MIXTURES 

Unlike simpler models such as the soft sphere mo del, the glass transition for the considered 
Lennard- Jones binary mixture cannot be described in the simple framework of well-known 
approximations of nonlinear integral equations for liquids: the Hypernetted Chain approxi- 
mation (HNC); the Mean Spherical Approximation (MSA); the Percus-Yevick Approximation 
(PY)il. Within all the classical approximation HNC, PY and MSA, the compressibility, com- 
puted as an integral of correlation functions: 



PXi = P + PP / d 3 rh(r) 



f33) 



diverges at a finite temperature Tql which depends on the chosen approximations [fig. 1] . We 
draw two conclusions from these computations. Our first conclusion is that this singularity 
corresponds to the liquid-gas phase separation which is typical of simple systemsEl. To con- 
vince the reader, in [fig. 2] we show the phase diagram in the plane p — T computed within 
HNC approximation, where the region of phase coexistence is shown. 

Hence, our claim is that HNC and PY, in this model, badly overestimate the real liquid- 
gas phase separation temperature, Tql, while MSA approximation allows to get far closer to 
this transition. This is confirmed by our numerical simulations (see below) placing Tql for 
p = 1.2 at ~ 0.3. 
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FIG. 1. Compressibility, %i as a function of temperature, T, computed by HNC (filled line), PY 
(dashed line) and MSA (dotted line) approximations. 
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FIG. 2. Projection of the phase diagram in the p — T plane computed by HNC approximation. 
In the region on the left of the line, a thermodynamically stable homogeneous phase does not exist 
and the system undergoes a liquid-gas phase separation. 

A simple procedure for improving the integral equations has been proposed by Zerah and 
HansenM It mixes the HNC and MSA closures by means of a single parameter, a, that is 
chosen in order to reduce thermodynamic inconsistencies between the two different routes for 
computing the compressibility: 



Briefly, we recall that HNC approximation consists of neglecting the 'bridge' diagrams 
in the cluster Mayer expansion, while MSA treats the attractive term of the potential as a 
correction to the repulsive term. 

The 'mixed' integral equation which we have resort to is, then: 




(34) 



where the pressure, P, is computed via the virial equation: 




(35) 



ee'p 2 q^(r) = exp (-(3V { %(r)) 1 + 



eMf e€ '(r){w~'(r)-(3V { %(r)}}-l 
f ee '(r) 



(36) 
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where V^)( r ) an d V(R)( r ) are > respectively, attractive and repulsive contribution to potential: 



_(V(r)-V K 
V(R) = { o 

V W = { V(r) 

and the function W(r) is h(r) — c(r), with c(r) the direct correlation function. 

In the Zerah and Hansen scheme the function f ee (r) is introduced, allowing continuously 
interpolation between HNC and MSA, simply by choosing: 

/"'(r) = l-exp( —). (38) 

Let us stress that for /(r) — > 1 the closure (^) reduce to HNC, while in the opposite limit, 
f(r) — > 0, MSA is recovered. The choice of (|38| ) implies that we obtain correlation functions 
which are HNC-like at large r and MSA-like at short r. 

We have chosen to utilize this scheme in a different fashion, since we are mainly interested 
to obtain pair correlation functions which very accurately reproduce the observed liquid be- 
haviour in order to study the effects of the used approximations in describing the liquid-glass 
transition. The HNC closure is exact at very high temperatures and, within this model, 
significant deviations from it of the numerical data of the energy appear yet in the region 
(3 G [10 _4 ,0.1]. On the other hand, the MSA approximation turns out to be more correct 
when the liquid-gas coexistence region is approached, both in reproducing the numerical data 
and in reducing the inconsistencies between the two different compressibilities. We have cho- 
sen as a starting value a = 1.5 at (3 — 0.1, which fits very well the numerical data in the 
whole region (3 G [10 -4 , 0.1] and we increased it linearly in such a way to have a MSA value 
(a = 5) at the liquid-gas transition. 

The result is a very good agreement between internal energy computed with the g(r) ob- 
tained by solving ( |36l) and the results of numerical simulations [fig. 3]. These will be discussed 
in detail in the following. Let us remind that the pure HNC and MSA closure have an error 
up to 10%, while, with interpolated integral equations, it is greatly reduced. 



■zL ' mm. 
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Hnc-Msa 
Simulations 



3.5 



FIG. 3. Internal energy vs. inverse temperature, j3. The points are obtained by simulated anneal- 
ing runs, the solid line by interpolated integral equation 



VI. RESULTS AND DISCUSSION 

We have studied a binary mixture where particles interact by means of a Lennard- Jones 
type potential. This model has been introduced and extensively studied as a good glass for- 
mer. For density p = 1.2, numerical simulations show evidence for a 'dynamical' transition 
at temperature Tp = 0. 4350. This kind of transition does not correspond to a real thermo- 
dynamic transition with singularities in thermodynamic quantities, but rather seems related 
to changes in the relaxation processes of the system in the free-energy landscape. Mean field 
theory of supercooled liquid, Mode Coupling TheoryS, describes this dynamical transition. 

Complementary to this 'dynamical' picture of the glass transition is the investigation of the 
free-energy of a 'molecular' liquid, where 'molecules' are formed by m replicas of the particles 
in the original system, as recently proposedlhi and studied in this paper for a Lennard- Jones 
model. 

Following this strategy, we show that at the temperature T K = 0.32 (see the figures below) 
the complexity, i.e. the entropy due to the large number of minima in potential energy as a 
function of particle positions, vanishes. In our 'thermodynamic' picture, Tk corresponds to 
the Kauzmann temperature, where the supercooled liquid undergoes a true thermodynamic 
transition and the system enters a new phase, where it behaves as an amorphous solid. 



A. Thermodynamic quantities 

From the thermodynamic point of view, the transition is second order. Indeed, the free 
energy and the internal energy have no singularity [fig. 4], while one can observe that the 
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specific heat [fig. 5] shows an evident 'jump' at Tk, remaining close to 3/2, in agreement with 
the Dulong-Petit law, throughout the glass phase. 




0.15 0.2 0.25 035 04 045 88 °- 15 02 025 -p 03 °- 35 04 045 

FIG. 4. The free energy (a) and the internal energy (b) of the Lennard-Jones model as a function 
of the temperature, in both the liquid and the glass phase. 
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FIG. 5. The specific heat of the Lennard-Jones as a function of the temperature. 

B. New thermodynamic parameters 

In the framework of replica theory of structural glasses two new parameters are introduced, 
analogous to those used to describe the glass phase in spin-glass models: the parameter m, 
plotted in [fig. 6a] and cages size, A + and A~ , plotted in [fig. 6b]. Both are nearly linear 
with temperature. In particular, this means that the effective temperature, f3 m, is always 
close to the transition value, so in our theoretical computation we need only the mean values 
of observables at temperatures T > T GL . 
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FIG. 6. m (a) and the cage radii A + and A (b) as functions of the temperature. 



C. Comparisons with MC simulations 

Before discussing the numerical results on the complexity and the evaluation of the tran- 
sition temperature, let us pay some attention to the possible presence of a liquid-gas phase 
separation at a non zero Tql- We identify this temperature as the one where the liquid 
pressure Pu q equals the gas pressure P gas - At low enough temperatures, the coexistence is 
between liquid and gas at very small density, where P gas ~ p/ (3 is negligible. Therefore, as a 
first approximation, we evaluate Tql as the temperature where the liquid pressure becomes 
compatible with zero. This pressure is computed via the virial equation: 

- = - " — fe>-«— (>«)V (39) 

We present in [Fig 7] data on the behaviour of ^(rdV/dr), that seems very well fitted by the 
power law aT 3 / 5 + b, usually encountered in liquidsS. By extrapolating data at temperatures 
lower than the ones where we succeeded in thermalizing the system, we get a definitely non 
zero T GL ~ 0.3. 
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" ? 0.2 0.4 _j_ 0.6 0.8 1 

FIG. 7. The numerical evaluation of Tql ~ 0.3 as the temperature where the extrapolation of 
data on ^(rdV/dr) equals 1//?, i.e. the liquid pressure becomes compatible with zero. 

Though more extensive investigations are needed to confirm this result, our data are in 
good agreement with the existence of a liquid-gas phase separation. Moreover, they strongly 
support our hypothesis that approximation schemes like HNC and PY overestimate the tran- 
sition temperature. Therefore, introducing the Zerah and Hansen scheme of interpolation is 
necessary to analytically evaluate liquid quantities near the transition point. Let us recall 
that this closure gives Tql — 0.25, a slightly lower value than numerical simulations. 

Details of the simulations will be discussed below but we would like to emphasize here 
that the model we considered is the LJ binary mixture described in (|j) without any cutoff on 
the potentials, taking into account the small correction due to finite size effects. The presence 
of a cutoff (it is usual to take V et (r) = for r > 2.5cr ++ ) would obviously influence both the 
energy and the pressure behaviour and we expect a lower value for Tq L . 

In order to evaluate the complexity and the thermodynamic temperature by simulations, 
we consider the numerical approach suggested by the harmonic resummation scheme, first 
utilized for the soft spheres binary mixture modeP. We measure both the liquid and the 
amorphous solid entropy, whose difference gives S. Tr is the temperature where the entropies 
become equal. 

The liquid entropy can be obtained by numerically integrating the energy density 

S Uq (P) = (3 (e H M - f liq (p)) = S? iq + (3 e H M - f d(3' e liq ({3') (40) 

•/ 

where Sf iq is the entropy of the perfect gas in the (3 —>■ oo limit, i.e in the binary mixture case 

Sf iq = 1 - logp - clogc - (1 - c) log(l - c). (41) 

As it has been discussed in the previous papeiH, the numerical computation of the 'har- 
monic solid' entropy is a more subtle task. One has 
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S sol ((3) = ^(1 + log(27r)) - ^ (Tr log(^M)) 



(42) 



but the measure of S so i from the 'instantaneous' Hessian is complicated by the negative eigen- 
values, always present beyond the mean field approximation, though their number decreases as 
exp(— 1/T) at low temperatures and it is expected to be negligible below the Mode Coupling 
temperature0'il. 

We consider three different ways of evaluating S so i: 

• We measure it from only the M pos positive eigenvalues 



q(a) 
°sol 



d 

2 



1 + log( 



,27T, 



Mr 



pos j = i 



(43) 



We consider also the absolute values of the negative ones (we are disregarding both here 
and in the previous case the very few |A| < 10~ 4 ), 



S. 



(b) 
sol 



1 + log( j) 



<1N 



,*vS'° g|A - 



(44) 



• Starting from an equilibrium configuration at a given T value, we perform a Monte 
Carlo run at T = 0, allowing only small displacements to each particle. The percentage 
of non-positive eigenvalues becomes very rapidly < 2% in the whole temperature range 
considered and, correspondingly, we get compatible results on when using only the 
positive eigenvalues or also the negative ones. 

As shown in [fig. 8], there is not a small difference between S^] and Sg% particularly at 
temperatures T ~ Tp. This is also related to the fact that nearly all the negative eigenvalues 
are less than one in absolute value. On the other hand, definitely appear nearer to the 
'instantaneous' Hessian entropy evaluated by taking all the eigenvalues in absolute value. 
This has also been observed! in the case of the soft sphere binary mixture and it suggests 
that is a more reasonable estimation of the solid entropy. 

The other problem that we face is the well-known and difficult task of thermalizing glass- 
forming liquids at low temperatures. We choose to perform a simulated annealing run of 
a rather large system, using data on the liquid energy down to the temperature where the 
equilibrium was still reachable in a reasonable CPU time (/3 ~ 2.5). 

Then we extrapolate the liquid entropy behaviour at lower temperatures by fitting data 
in the interval (3 G [1, 2.5] with the power law 

S liq (T) = aT- 2 " > + b. (45) 

This corresponds to the assumption of a low temperature liquid energy behaviour e ~ T 3 / 5 , 
a theoretical prediction that seems to be well-followed in many simple liquidsll and that we 
find in very good agreement with our data, both in this case and in the previously discussed 
data on the pressure. 
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In more detail, we performed a simulated annealing run of a system of N — 260 particles, 
in a cubic box with periodic boundary conditions, starting from p 1 ^ = 0.02 and performing 
up to 2 22 MC steps at each A/9 1//4 = 0.02, the maximum shift 5 max allowed to each particle in 
one step being chosen in order to get an acceptance a ~ 0.5. The energy and its fluctuation 
were measured in the last half of each run at a given (3 value while the pressure was obtained 
from about 200 configurations in the same interval. We take into account the small finite size 
effects by adding to the energy and pressure the appropriate corrections evaluated analytically 
by means of J A v (r)d 3 r, where A is the volume obtained as difference between the volume of 
the whole space and the volume of the box containing the particles. Just for decreasing 
the error on the evaluation of Su q , we fit the very high temperature data on the energy, up 
to (3^ 4 = 0.1, by using (3 3 / 4 e((3) = a(3 1 / 2 + bf3 1 / 4 + c, obtaining correspondingly f([3 ) = 
4a/5g / ' 4 /3 + 2b/3^ 2 + 4c/3o, that turns out to be perfectly compatible with the analytical value 
(i.e. we are still in the region where no differences are observable between numerical data 
and analytical approximations). The integration is subsequently performed by interpolating 
the simulation data with a standard numerical subroutine in order to get a result which is 
independent of the integration interval. 

In order to evaluate S^j and Sf^ we considered 16 different configurations in the last half 
of the run at each /3-value, while S^j was measured from the configurations obtained from 
these by 10000 MC steps at T = (starting from 5 max = 0.05 and decreasing it up to 0.005 
during the run). 




1 1.5 2 2.5 3 3.5 

P 

FIG. 8. The entropy of the liquid (+) and the different evaluations (see text) of the amorphous 
solid entropy, (x), (*) and S^) (□), as functions of 0. The lines are the best fit to the 
power law, Su q = a(3 2 ^ 5 + b. 

In [fig. 8] we plot both Su q {(3) and the different measurements of S so i((3) considered. Sg b Q \ 
and are very close and give similar estimations of Tk ~ 0.3, whereas gives the slightly 
lower value T K ~ 0.26. It seems therefore reasonable to get the mean value T K ~ 0.28 as our 
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numerical evaluation of the Kauzmann temperature. 




FIG. 9. Analytical Su q (upper line) compared with the numerical values (+ points), and analytical 
S so i computed in the harmonic resummation scheme (lower line) compared with the three different 
numerical results, i.e. S$ (x)> ^2 (*), (□). 



We note that this value of the Kauzmann temperature, to be compared with the analytical 
Tk — 0.32, is close to the previously discussed numerical estimation of Tql- Our numerical 
data are inadequate for clarifying this situation, which is intriguing and deserves a more 
careful study. Here we only suggest that there could be a glass-gas phase coexistence in the 
low temperature region. 

Finally, [fig. 9] illustrates the comparisons between Su q and S so i as functions of the temper- 
ature T computed analytically (lines) and the same quantities determined by our numerical 
simulations (points), as discussed above. In this plot we added both the correct normaliza- 
tion constant, giving S = when the volume of phase space is equal to h 3 , and the thermal 
contribution which had been neglected in the previous figure. Obviously these terms do not 
affect the thermodynamical glass transition. We used the physical parameters appropriate to 
argon (<r ++ = 3.405 A, e ++ /kh = 125. 2K), plotting the entropies per particle. 
The reliability of the approximations introduced is clarified by these comparisons: 

• the numerical points of Su q are very close to the corresponding analytical quantity 
computed by means of the interpolation between HNC and MSA closures 

• the three different ways to numerically evaluate S so i are in various agreement with 
the analytical curve, but it is important to note that in the low temperature region 
(below the dynamical temperature), none of them is further than ~ 10 — 20% from the 
theoretical prediction. 
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